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Effects of impurity scattering on tunneling conductance in dirty normal- 
metal/insulator/superconductor junctions are studied based on the Kubo formula and 
the recursive Green function method. The zero-bias conductance peak (ZBCP) is a 
consequence of the unconventional pairing symmetry in superconductors. The impurity 
scattering in normal metals suppresses the amplitude of the ZBCP. The degree of the 
suppression agrees well with results of the quasiclassical Green function theory. When 
superconductors have d + is-wave pairing symmetry, the time-reversal symmetry is bro- 
ken in superconductors and the ZBCP splits into two peaks. The random impurity 
scattering reduces the height of the two splitting peaks. The position of the splitting 
peaks, however, almost remains unchanged even in the presence of the strong impurity 
scattering. Thus the two splitting peaks never merge into a single ZBCP. 

KEYWORDS: Andreev reflection, d-wave superconductor, zero-energy states, zero-bias conductance peal<, 
disorder, d + is-wave state 



§1. Introduction 

Charge transport in unconventional superconductors shows remarkable properties which can- 
not be expected in conventional s-wave superconductors. One of the striking effects is the zero-bias 
conductance peak (ZBCP) in normal- metal/unconventional superconductor junctions.^) In hybrid 
structures consisting of high-Tc superconductors, a number of experiments observed the ZBCP.^"^'') 
The zero-energy state (ZES)^^^ formed at a surface of the unconventional superconductor is respon- 
sible for the ZBCP.^^ The interference effect of a quasiparticle in the presence of the sign-change in 
the pair potential is the origin of the ZES.^^"^®-* The ZES is also responsible for the low-temperature 
anomaly of the Josephson current in unconventional superconductor junctions^^"^^^ and anomalies 
of the charge transport in ferromagnet/unconventional superconductor junctions. ^^"'^^^ 

Since the ZES is a result of the interference effect of a quasiparticle, it is sensitive to the 
time-reversal symmetry (TRS) of systems. Actually it has been confirmed in normal-metal/d-wave 
superconductor junctions that the ZBCP splits into two peaks by applying magnetic fields. ^^"^^-^ 
It has been also pointed out that the ZBCP splits into two peaks when the TRS is broken in 
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superconductors.^^' Theories showed a possibihty of d + is-^^^ or d + id'^'^^ -wave broken time- 
reversal symmetry states (BTRSS) at surfaces of superconductors. If such BTRSS's are stable at 
the surface of superconductors, they may bring significant influences on transport properties. ^^'^^^ 
Experimental results, however, are still controversial. Some experiments reported the split of the 
ZBCP in the absence of magnetic field, ^^"^^^ others did not show the zero-field splitting even in low 
temperatures.^' Moreover a experiment of the Josephson junctions concluded the absence 
of the BTRSS at interfaces of d-wave junctions. Thus it is very important to theoretically address 
the origin of the splitting. There are several factors which determine the magnitude of splitting of 
the ZBCP: the transmission probability of the junctions,^^) the roughness of the interface,^^) the 
impurity scattering in superconductors^^) and the impurity scattering in normal metals. Among 
them, effects of the impurity scattering in normal metals has not been studied yet. We note that 
effects of random potentials on the split of the ZBCP are also unclear within theoretical studies at 
present. Some theories based on the quasiclassical Green function method^^"^^) did not show the 
splitting,^^"^^) whereas others concluded the split of the ZBCP by a numerical simulation^^^ and 
an analytic calculation.^^) When the BTRSS is formed at the NS interface, a group of experiment 
predicted that the two splitting peaks may merge into a single ZBCP by the impurity scattering.^^) 

In dirty normal-metal/insulator /super conductor (DN/I/S) junctions, there are many theoret- 
ical studies (see a review^^)) on conventional s-wave superconductor junctions such as numeri- 
cal methods,^^^^) the random matrix theory^^' and the quasiclassical Green function ap- 
proach.^^"^^) These studies show a wide variety in line shapes of the conductance depending on 
the Thouless energy and the resistance of DN. In particular, the ZBCP is an important conclusion 
in the low transparent junctions even in the s-wave junctions. We note that the interference of 
a quasiparticle while traveling the diffusive metals causes the ZBCP in the s-wave junctions. Re- 
cently, the full resistance (R) of the disordered junctions is derived from a microscopic theory. 
By applying this theory to d-wave junctions, we show various transport properties which reflect the 
unconventional pairing symmetry in superconductors. When the a axis of high-Tc superconductors 
is oriented 45 degrees from the interface normal, R at the zero temperature can be expressed as 
Rd + RcJ^^ where Rd is the resistance of a diffusive metal and Rc is the resistance of a clean 
junction. The resistance of the disordered junction is given by the simple summation of the two 
resistance, which indicates the absence of the proximity effect. ^^^-"^^ At the present stage, however, 
there is no corresponding studies in DN/I/d-wave superconductor junctions with the BTRSS. 

In this paper, effects of the randomness on the ZBCP in DN/I/S junctions are studied numer- 
ically by using the recursive Green function method. We mainly consider the d + is-wscve pairing 
symmetry in superconductors. For comparison, we also discuss the conductance in the s- and d- 
wave junctions which are the two limitting cases of the d + zs-wave junctions. First we confirm the 
proximity effect in s-wave junctions. The total resistance at the zero-bias in disordered s-wave junc- 
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tions shows the reentrant behavior as a function of Rd, which is consistent with the quasiclassical 
theory. We also confirm that the resistance in disordered d-wave junctions increases monotonically 
with the increase of Ro^'^'^^^ since there is no proximity effect. The calculated results are consis- 
tent with the quasiclassical theory.^^^ In d + is-wave junctions, we show that the proximity effect 
appears when s-wave component becomes dominant. In d + is junctions, the ZBCP splits into two 
peaks because of the BTRSS and the impurity scattering reduces the height of the splitting peaks. 
The two splitting peaks, however, do not merge into a single peak. In the light of our theory, the 
single sharp ZBCP observed in experiments^' ^' is a strong evidence for the formation of the pure 
d-wave symmetry state at the interface. 

This paper is organized as follows. In Sec. II, the model and method are presented. The 
numerical results are shown in Sec. III. In Sec. IV, we summarize this paper. 

§2. Model and Method 

Let us consider a normal- metal / superconductor junction on the two-dimensional tight-binding 
model as shown in Fig. 1, where r = {l,m) labels a lattice site. We assume the periodic boundary 
condition in the y direction and the width of the junction is Mao, where uq is the lattice constant. 
The BCS Hamiltonian of the system is expressed as, 

Hbcs = -t [4,aV,a + H.c. 

{r,r'),a- 
r,a 

-EE[^*(^»^<TCr-a+^^-C-], (1) 
r r' 

where ct^a (cr.o-) is the creation (annihilation) operator of an electron at (/, m) with spin cj(=| or 
I) and /U is the Fermi energy. The summation ^(^^^/^ runs over nearest neighbor sites and t is the 
nearest neighbor hopping-integral. We assume that the on-site potential Vr is zero far from the 
interface in the normal metal and in the superconductor, (i.e., vi^^ = for Z < and L -|- 2 < Z). In 
the disordered region, 1 < I < L, the potential is given by random number uniformly distributed 
in the range of 

W W 

-— < vi,„. < —. (2) 

The insulating barrier is described by the potential VL+i,m = ^ins for all m. When a superconductor 
has the s-wave pairing symmetry, the pair potential is given by 

^ ' 1 : otherwise ^ ^ 
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When the a axis of a high-Tc superconductor is oriented by 45 degrees from the interface normal, 
the pair potential is given by 



A\r',r) = < 



Ad 




r lb 1, m = m' lb 1 
I' ziz l,m = m' ^ 1 
otherwise 



(4) 



We note that the tight-binding lattices do not correspond to the two-dimensional Cu02 plane in 
high-Tc superconductors. The tight-binding model represents the two-dimensional space. In our 
model, we introduce the pair potential between the next nearest neighbor sites to describe junctions 
under consideration. An alternative way is keeping the pair potential between the nearest neighbor 
sites and rotating the square lattice by 45 degrees. There are no essential differences between 
results in the two models when we focus on the formation of the ZES. This is because the ZES is 
a consequence of the d-wave symmetry of the pair potential. The BTRSS is characterized by the 
d + zs-wave symmetry and the pair potential is described by 



A'^+*'*(r',r) = A'^(r',r)cosa-biA'*(r',r)sina, 



(5) 



where 0<Q;<7r/2isa parameter which characterizes the degree of the time-reversal symmetry 
breaking. When s-wave component appears only at a surface of a superconductor, we describe the 
pair potential as 



A'i+''{r',r) =A'^(r',r)tanh ( | 



+ iA^(r',r)|l-tanh(^|j|, 



(6) 



where ^ is comparable to the coherence length of d-wave superconductors.^^' In a normal metal, 
the pair potential is taken to be zero. 

The BCS Hamiltonian in Eq. (1) can be diagonalized by applying the Bogoliubov transforma- 
tion. 



J 

-r,l 



E 



(7) 



Mr) Kir) 

where "ft^a {1u,(t) is creation (annihilation) operator of a Bogoliubov quasiparticle. We omit the 
spin index of the wavefunction (n^,Wjy) since we do not consider the spin-dependent potential. In 
this way, the Bogoliubov-de Gennes (BdG) equation is derived on the tight-binding lattice. We 
solve the BdG equation by using the recursive Green function method. ^^'^^-^ The Green function 
is defined in a 2M x 2M matrix form 



G(M') = E 



(8) 
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where E^, is the eigen value of the BdG equation, ^ui^) is the vector with 2M components and 
m-th (m + M-th) component is Ui,{l, m) {vy{l, m)). The differential conductance of the junction is 
calculated to be 



GNsieV) =-t^Tr' 6(1,1 + 1)0(1,1 + 1) 
n I 

+ G(l + 1, l)G(l + 1,1)- G(l, l)G(l + !,/ + !) 



(9) 



with E = eV, where V is the bais voltage applied to the NS junction. When j4 is a 2M x 2M 
matrix, Tr'[4] means the summation of X]m=i ^m,m- The conductance in this method is identical 
to that in the well known conductance formula.®^' By averaging over different configuration of 
randomness, the mean value of conductance < Gns > is obtained. The conductance in normal 
conductors G^ can be also calculated in the conventional recursive Green function method. '^^^ 
The advantage of the method is wide applicability to various systems such as, the clean (ballistic) 
junctions, the dirty (diffusive) junctions and the junctions in the localization regime. ^^-^ 



Insulator 



Normal metal 
Clean region Disorderd region 1 



Superconductor 




m=M 



m=\ 



L L+lL+2 



Fig. 1. The normal-metal / superconductor junction is illustrated, where / and m are the lattice indices in the x and 
the y direction, respectively. The width of the junction is M. We introduce the disordered region with length L 
in a normal metal near the interface ( gray circles) and the insulating barrier at Z = L -I- 1 ( filled circles). 



§3. Results 

Throughout this paper, we fix the Fermi energy at = —l.Ot, the amplitude of the pair 
potential in the d-wave component at A^^ = O.OOlt and the impurity potential at W = 2.0t. In 
order to determine the mean free path at the Fermi energy, we first calculate the normal conductance 
as a function of the length of the disordered region L in Fig. 2(b). The ensemble average is carried 
out over 500 samples with different random impurity configurations. When L is small, (Gm) ~ 
(2e^//i)M. Thus the disordered region is in the quasiballistic transport regime and {Giy)L/M 
is proportional to L as shown in L < 50. In the diffusive regime, {Gn)L/M coincides with the 
conductivity which is independent of sample size. Thus 60 < L < 150 in Fig 2(b) corresponds 
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to the diffusive regime and the mean free path is estimated to be about 6 uq. For large L, the 
disordered region is in the locahzation transport regime and the conductance decreases with L in 
proportional to exp{— L / $,al) i where ^al is the localization length of an electron. 



(a) 

Normal metal 
clean disorderd clean 




Fig. 2. In (a), the system used in a simulation is schematically depicted. The normal conductance is plotted as a 
function of the length of the disordered region in (b), where W = 2.0t, \i — — l.Ot and M = 32. When the 
disordered region is in the diffusive regime, the vertical axis in (b) corresponds to the conductivity. 

In Fig. 3, we show the resistance at the zero-bias in DN/I/S junctions as a function of L. We 
assume d + is pairing symmetry in superconductors as shown in Eq. (5), where Ag = A^. The 
results for a = 7r/2 correspond to the resistance in s-wave junctions. The resistance for a = 7r/2 
first decreases with the increase of L. Cooper pairs penetrate into normal metals and have a finite 
amplitude. This is a source of the proximity effect. As a consequence, the reflectionless tunneling 
suppresses the resistance at the zero-bias. For large L, the resistance increases linearly with L. 
The reentrant behavior of the resistance for a = 7r/2 is the direct consequence of the proximity 
effect. The results for a = correspond to the resistance in d-wave functions. The resistance does 
not show the reentrant behavior and increases almost linearly with the increase of L. In this case, 
the proximity effect is absent. The amplitude of the Cooper pairs becomes almost zero because 
the contribution of the pair potential with positive sign and that with negative sign cancel out 
with each other. For a = O.lvr, the proximity effect appears because of the s-wave component in 
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the pair potential and the resistance show the weak reentrant behavior as a function of L. For 
a = O.Svr, the reentrant behavior becomes rather remarkable and the results are close to those in 
the s-wave junctions. According to the quasiclassical Green function theory, the resistance of the 
cZ-wave junctions can be simply given by 



where Rd is the resistance in DN and Rq is the resistance of the clean junction. The summation of 
these two resistance gives the full resistance of the disordered junctions because there is no proximity 
effect. The absence of the proximity effect is responsible for the anomaly in the Josephson current 
through diffusive metals. ^'^'^^^ 



Fig. 3. The resistance of the DN/I/S junctions is shown as a function of the resistance in DN (Rd), where Vins = 5t 
and M = 32. The pair potential in superconductors is given in Eq. (5). 

In Fig. 4, we show the conductance in clean N/I/S junctions as a function of the applied bias 
voltage, where Vins = 5t and = A^- For a = O.Svr, junctions become s-wave junctions and the 
conductance has two peaks at E = ibAg which reflect the singularity of the density of states in 
s-wave superconductors. Junctions with a = correspond to d-wave junctions and the results show 
the ZBCP owing to the ZES at the junction interface. For a = O.Olvr, the ZBCP splits into two 
peaks because the TRS is broken by the s-wave component in the pair potential. ^^-^ The degree of 
the splitting increases with increasing a. The calculated results basically remain unchanged when 
we assume the s-wave component only at the interface. Since the splitting was observed in several 
experiments, ^^"^^^ the zero-field splitting of the ZBCP has been believed to be an evidence of the 
BTRSS. On the other hand, other experiment^' ^'^'^'^^'^'^^ did not show the zero- field splitting. Thus 
the experimental results are still controversial at present. There are some reasons which explain 
the contradiction in experiments. One of them is the sample quality of the junction. Actually a 
group of experiment insists that the zero-field split would be found when the sample quality is good 
enough and the measurement is done in sufficiently low temperatures. If this prediction is true, the 
impurity scattering might unify the splitting peaks. We try to check this prediction as follows. 



R = Rd + Rc 



(10) 



200^ 



a 




-A- 
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Fig. 4. The conductance in clean normal-metal / superconductor junction is plotted as a function of the bias voltage, 
where Vins = 5t. The s- and d-wave junctions are correspond to a = 7r/2 and 0, respectively. 

We introduce the impurity potential in normal metals of d + is-wave junctions. In Fig. 5, the 
conductance is plotted as a function of the bias voltage for several L, where Vins/t = 2, 7 and 10 
in (a), (b) and (c), respectively. The s-wave component is introduced only at the NS interface and 
^ in Eq. (6) is chosen to be 10 lattice constants. The amplitude of the s-wave pair potential is 
fixed at = 0.2Aii- When the transparency of the interface is high, the ZBCP in clean junctions 
does not split into two peaks even in the presence of the s-wave component as shown in (a). The 
splitting can be seen in the disordered junctions as shown in the results for finite L. Thus the 
impurity scattering in normal metal seems to assist the splitting of the ZBCP. At the same time, 
the amplitude of the conductance decreases with increasing L. In highly transparent junctions, 
the calculated results indicate the opposite conclusion to the experimental prediction. When the 
transparency of the junction becomes low, the ZBCP splits into two even in clean junctions as show 
in (b) and (c). The amplitude of the conductance decreases with the increase of the resistance of 
normal metals as well as those in (a). It is important that the degree of the splitting are not 
changed by the impurity scattering. In low transparent junctions, our results are also contradict to 
the experimental prediction. We conclude that the impurity scattering in low transparent junctions 
does not merge the splitting peaks into a single ZBCP. 

§4. Conclusions 

We numerically calculated the differential conductance in dirty normal metal/ insulator/ su- 
perconductor junctions by using the recursive Green function method. In superconductors, we 
assumed d + is pairing symmetry which breaks the time-reversal symmetry of the junctions. In 
pure s-wave (pure d-wave) junctions, we confirmed the presence (absence) of the proximity effect, 
which agrees with the quasiclassical Green function theory. In the case of the d + is symmetry, 
the proximity effect recovers when the pair amplitude of the s-wave symmetry becomes dominant. 
In highly transparent clean junctions, we found the single ZBCP even in the absence of the time- 
reversal symmetry. The ZBCP, however, splits into two peaks owing to the impurity scattering. 
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. 5. The conductance in DN/I/S is plotted as a function of the bias voltage, where M = 32, ^ = 10 and As = 0.2A, 
The potential of an insulating layer are chosen to be Vins/t = 2, 7 and 10 in (a), (b) and (c), respectively. 
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The s-wave component splits the zero-bias conductance peak into two peaks in low transparent 
junctions even in the clean limit. The amplitudes of the splitting peaks decreases with the increase 
of the impurity scattering in normal metals. The peak position, however, remains unchanged even 
in the presence of the strong impurity scattering. Thus we conclude that the impurity scattering 
does not merge the two splitting peaks. 
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